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Abstract 



The large- limit of the expectation values of the Wilson loops in the fun- 
damental representation, corresponding to two-dimensional U(A) Yang- 
Mills and generalized Yang-Mills theories on a sphere are studied. The 
' behavior of the expectation values of the Wilson loops both near the 

, critical area and for large areas are investigated. It is shown that the ex- 

' pectation values of the Wilson loops at large areas behave exponentially 

l/^ ' with respect to the area of the smaller region the boundary of which is 

. the loop; and for the so called typical theories, the expectation values of 

the Wilson loops exhibit a discontinuity in their second derivative (with 
respect to the area) at the critical area. 

^ ; 1 Introduction 

Two-dimensional Yang-Mills theory (YM2) and generalized Yang-Mills theories 
^ . (gYM2's) have been a subject of extensive study during recent years [1-17]. 

' These are important integrable models which can shed light on some basic fea- 

tures of QCD4. Also, there exists an equivalence between YM2 and the string 
theory. It was shown that the coefficients of the 1/A^ expansion of the par- 
tition function of SU(iV) YM2 are determined by a sum over maps from a 
two-dimensional surface onto the two-dimensional target space. 

The partition function and the expectation values of the Wilson loops of YM2 
have been calculated in lattice- [1,18] and continuum-formulations [4,12,13]. 
The partition function and the expectation values of the Wilson loops of gYM2 's 
have been calculated in [10, 11]. All of these quantities are described as sum- 
mations over the irreducible representations of the corresponding gauge group. 
In general, it is not possible to perform these summations explicitly. For large 
gauge groups, however, these summations may be dominated by some specific 
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representations, and it can be possible to perform the summations explicitly. 
There are other physical reasons as well (for example the relation between large- 
N YM2 and the string theory), that make the study of the YM2 and gYM2's 
for large groups important. 

In [19], the large- iV limit of the U(iV) YM2 on a sphere was studied. There 
it was shown that the above mentioned summation is replaced by a (functional) 
integration over the continuous parameters of the Young tableaux corresponding 
to the representation. Then the saddle-point approximation singles out a so- 
called classical representation, which dominates the integration. In this way, it 
was shown that the free energy of the U(7V) YM2 on a sphere with the surface 
area A < Ac = has a logarithmic behavior [19]. In [14], the free energy was 
calculated for areas A> ir'^, from which it was shown that the YM2 on a sphere 
has a third-order phase transition at the critical area Ac = tt^, like the well 
known Gross- Witten-Wadia phase transition for the lattice two dimensional 
multicolour gauge theory [20,21]. The phase structure of the large- A'' YM2, 
generalized YM2's, and nonlocal YM2 on a sphere were further discussed in 
[15,17,22,23]. It is also seen that for surfaces with no boundaries, only the 
genus-zero surface (i.e. the sphere) has a nontrivial saddle-point approximation, 
and all other surfaces have trivial large- behavior [19]. 

In [24-26], the expectation value of the Wilson loops corresponding to the 
large- A/' U(A'')-YM2 on a sphere was investigated. The study was again based 
on the fact that in the large- TV limit, a certain representation of the group is 
singled out in the partition function and the expectation value of any quantity 
is essentially the value of that quantity calculated at that representation. Using 
this, the behavior of the expectation values of the Wilson loops for large- A/' YM2 
was investigated for small and large areas. 

In this paper we want to study the expectation-value of the Wilson loops 
in the fundamental representation, corresponding to large- A' YM2 and gYM2's 
on a sphere. In section 2, the expectation values of a Wilson loop is obtained 
in terms of the the area enclosed by the loop and the density of the dominant 
representation. In section 3, the behaviors of the free energy and the expectation 
values of the Wilson loops for areas near the critical value are investigated by 
some general arguments. In section 4, we justify explicitly the results of section 
3 for two special cases YM2 and G{z) = gYM2 models. Finally, in section 5 
the large-area behavior of the expectation value of the Wilson loops is studied. 

2 The expectation values of the Wilson loops 

Following [11, 17], a generalized U(A'') Yang-Mills theory on a surface is charac- 
terized by a function A: 



p 




(1) 
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where R denotes a representation of the group U(7V), a^'s are constants, and 
Cfc's are the Cashnirs of the group defined through 



N 



Ck^J2i(n, + N-if-{N-i)% (2) 



rii's are nonincreasing integers characterizing the representation. It is assumed 
that p is even and ap > 0. For simplicity, from now on it is further assumed 
that aU Ofc's with odd fc's vanish. The partition function of such a theory on an 
orientable genus g surface with n boundaries is 



{A; C/i, ...,[/„) = ^ 4-'^-" XRiUi) ■ ■ ■ XR{Un) exp[-A A(i?)], (3) 



where UiS are the holonomies on the boundaries, A is the area of the surface, 
da is the dimension of the representation i?, and XR is the character of the 
representation R. 

Consider a simple (not self-crossing) loop, which divides a sphere into two 
regions of areas Ai and A2. The expectation value of the Wilson loop corre- 
sponding to this loop is (following [11,17,24]) 

^r^^^^^^)= Z^{aI+A2) J ^oAAi,U)ZoAA2,U-')j^^Xr{U), (4) 

where r is a representation of U(A^), Zo{A) is the partition function on the 
sphere, and the integration is done over the group. Eq. can be rewritten as 

Wr{Ai,^2) = , y E ^« '^s {R, r\S) exp[-Ai A(i?) - A, A(5)], 

(5) 

where the summations run over irreducible representations ofV{N), and (i?, r\S) 
is the number of representations S in the tensor product of the representations 
R and r. From now on, the representation r is taken to be the fundamental 
representation, the dimension of which is N. 

For large N, the summation on R is dominated by the so called classical 
representation, which maximizes the product (Ir exp[—Ai A(i?)], as was shown 
in, for example, [17]. To obtain this representation, it is convenient to introduce 
the new parameters 

i 

h{x) := —n{x) — 1 -|- x, (6) 

the density 

pW :^ , (T) 
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and the function 



G(z) j^a.l-z)^ 



(8) 



Then the density corresponding to the classical representation is characterized 

by 

A f 

— G{z) — dy p{y) In \z — y\ ~ const., iff p{z) ^ 1 and — a < z < 

dz p{z) = 1, (9) 



where a is a positive number to be determined through the above conditions. 
Then, using arguments similar to those used in [24], one arrives at 



— sin[7rp(z)] exp -P / 

-a ^ V J-a 



dy^ + A2 G'{z) 
z-y 



where P means Cauchy's principal value, or 



f dz 

W{Ai,A2) = - f TT- e^V 
Jc 27ri 



dy-^+^2 G\z) 

z-y 



(10) 



(11) 



G'{z) denotes the derivative of G{z), and C is a counterclockwise contour en- 
circling the real segment [—a, a]. Defining the function H of a complex variable 
z as 

p(y) 



H{z) 



dy 



z-y 



(as it was defined in [17]), one can rewrite as 

W{A^,A2) = - £ 1^ exp [-H{z) + A^ G'{z)] , 



(12) 



(13) 



It was seen in [17] that H is analytic except on the real segment [—a, a]. In 
Yang-Mills case where G{z) = (l/2)z^, leads to one derived in [24,25]. 



3 Free energy and the expectation values of the 
Wilson loops near the critical area 

As it was discussed in [14, 17], if the area of the sphere is less than the critical 
value then the density p is always less than 1 and the free energy of the 
system is a smooth function of the area. For areas larger than Ac, however, p 
becomes equal to 1 on certain segments. This produces a discontinuous behavior 
in the area-dependence of the free energy, at A = A^. There was an argument 
in [23] to obtain the degree of this nonsmoothness , i.e. the order of this phase 
transition. However, the argument is problematic and, as will be shown here, 
leads to incorrect results at some cases. 
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Let us denote the correct density for A > (the strong region) by ps 
and the density satisfying the saddle-point equation, but not subject to not 
exceeding 1, by pw Define and Ls as the regions where ps is less than one 
or equal to one, respectively; and L'^ and L'^ as the regions where pw is less 
than one or greater than one, respectively. The length of Lg is denoted by 26, 
so the length of L'^ is of the order 2b. We have 

^G'{z) = P fdyfM, ,eL^, (14) 
2 J z-y 

4g'(z) = P /dy^, zeL'^UL'^. (15) 
2 J z-y 

where ps {pw) has been defined zero outside U Lg [L'^ U L'^,). Defining 

a :— niax(pu,) — 1, (16) 

and using the ansatz 

Pso{y) ■■= mm[p^{y), 1] (17) 

as an approximate solution to the ea. (|14|l . it is seen that the difference of the 
left-hand side with the right-hand side calculated at ps = Pso, is of the order a 
when the distance of z from Lg is of the order 6, and of the order a b when the 
distance of z from Lg is large compared to b. To see this, one notes that 

Di^) ■.^^G'{z)-F UyP^'^y^ 



2 ' ' J z-y' 

J z-y 
=P / dy PjM^^ , e L'^, (18) 

Jl'^ z-y 

The numerator of the integrand is of the order a, and the length of the inte- 
gration region is of the order b. If the distance of z from Lg (or L'J is of the 
order b, then the denominator of the integrand is of the order b, and the integral 
would be of the order a. If the distance of z from Lg (or L'g) is large compared 
to 6, then the denominator of the integrand is of the order one, and the integral 
would be of the order a b. 
Defining Spg through 

pg =: pgo + Spg, (19) 

it is seen that 

Spsiy) 



D{z)=V Ay^^^, zeLu,. (20) 
J z-y 

Sps is at most of the order a. But 6p cannot be of order a everywhere. In fact, 
if it is of the order a in a region large compared to b, then the right-hand side 
would be of the order a for some points z the distance of them from Lg (or Lg) 
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is large compared to b. But D(z) is of the order a b when the distance of z from 
Ls (or L'g) is large. So, 

where L" is a region around Lg, the length of which is of the order b. 
Using lEU and 

:= / dy (22) 

J z-y 

it is seen that 

H,{z) - H^z) ^ f dy l^P^ + f dy = o(a), \z\ » a, (23) 

Jl, z-y Jr\l, z-y 

where a is of the order of the length of the region in which the integrand does 
not vanish. 

For large values of z, the difference {z)—Hw{z) is an analytic even function 
of b [17]. This means that this difference can be expanded like 

H,{z) ~ H^{z) = Cfe(z) b""" + cu+i{z) b^'+'^ + ■■■ , (24) 

where fc is a positive integer. If at the point attains its maximum, the 2m'th 
derivative of pw is its first nonvanishing derivative, then 6^™ is of the order a. 
As the difference Hs{z) — H^{z) is o{a), it turns out that k must be greater 
than m. So one arrives at 

Hs{z) - H.^{z) - 62m+2^ foj, large z (25) 

or 

Hs{z) - H^{z) - for large z. (26) 

The free energy F is defined through 

^:=-]^ln^, (27) 

and one has [17] 

F'{A)^ f Ayp{y)G{y). (28) 
Using (|12|1 . one can rewrite this as 

^^'(^)=£ ^i?(z)G(z), (29) 

where Coo is a large counterclockwise circle. Now, using (|26|l . it is seen that 

F^{A) - F^,{A) - (A ~ ^^)(2;-i)[i+(i/™)]^ (30) 
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and from that 

Fs{A) - F^{A) - (A - (31) 

where the {21 — l)'th derivative of a with respect to A is the first nonvanishing 
derivative of a with respect to A, so a{A) ^ {A — Ac)^^^^. It is seen that in 
general the last equation differs that found in [23]. The two, however, coincide 
if / = 1. The flaw in the argument of [23] is that from the fact that the integral 
of the difference ps — on is of order a b, it was deduced that ps — Pw 
itself is of order ab, which is not true, as seen from H21|) . However, for typical 
theories where Z = m = 1, it is still true that the system exhibits a third order 
phase-transition at A — Ac- 

Using H13|) and H26|) , one can also investigate the nonsmoothness behavior of 
the expectation values of the Wilson loops. Using (|f 3|l . It is seen that 

Ws-W^^iA- ^^)(2;-i)[i+(i/™)]^ (32) 

so that for typical theories, there is a second order discontinuity in the ex- 
pectation values of the Wilson loops. For example for all gYM2 theories with 
G{z) — flfe z'^''{k > 1), including YM2, it has been shown that [23] 

a{A) = i^-J 1. (33) 

This shows that the derivative of a with respect to A never vanishes, or Z = I. 
Also the second derivative of pw is negative at the maximum of p^ , which gives 
m = 1. So for all theories with G{z) = at z^^ , 

Ws - r^{A- A,f . (34) 



4 Explicit evaluation of the order of discontinu- 
ity of Wilson loops for YM2 and the gYM2 
with G{z) = 



From H12(l , and using the fact that p is even if G is even, one has 



, ^ 1 1 
Hiz) = - + - 



Using this and H13|l . one arrives at 
f dz 

Ws-W^^ i — exp [-H^{z) + A2 G'{z)] [Hs{z) - Hjz) - 

"A(y2) ^ A(y4) 



= exp[-H^(z)-t-A2G'(z) 



(35) 



(36) 
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where 



A(y") :=(2/"). - (2/")», 

dyy^p{y)[ps{y)-pM- (37) 



It is seen that Wg — Ww is proportional to A(?/^). For small values of A2, 
one can expand exp[AG"(z)] and it is seen that only a finite number of terms 
in are needed to obtain the behavior of Ws ~ Ww (In the product of 
[G'(2;)]" \ils(z) — iJ^(z)], one needs to keep only terms of for fc > —1.) 

For example, if G{z) = (1/2) z^, using the first two terms one can obtain the 
behavior of Wg — up to order For this theory, one has (using [14]) 



2 /A-Ar 



r2 



Ac 



(38) 



2 



and 



which shows that leading term of Wg — W^ is proportional to {A — Ac)'^, at least 
up to order (^2)'*- 

As a second example, one can take G(z) = z"^. Following [17], one can 
expand Hg — near the critical area, and expand it in terms of the powers of 
from which A(y^) (the coefficient of z~^) is calculated to be 

551 fA-Ar'^'^ 



which shows that the leading term of Wg — is proportional to {A — AcY , at 
least up to first order in A2. 



5 The expectation values of Wilson loops for 
large areas 

To study the Wilson loop on the plane {A — > 00,), one notes that for large values 
of A, the first equation of the set Q cannot be fulfilled with finite p. This shows 
that as A tends to 00, the density p becomes equal to one everywhere in (—a, a). 
Then the second equation shows that a tends to (1/2) as A tends to 00: 

lim p{z) — 1, —a < z < a, 

A^oo 

lim a = i. (41) 

A—>oo 2 
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This shows that 



2+1 

lim H{z) = In y- 



A^oo 



(42) 



For A large but not infinite, one can use H12|) to expand H{z) around the value 
of H for infinite A. The result would be 



E 



{z - z,Y 



(43) 



where n runs over positive integers and z^'s are points in [—\, \\ around them 
p differs from one (for large but not infinite areas) . Among these points are — ^ 
and ai^„'s should of course tend to zero as A tends to infinity. Using the 
above expansion in ()13() . one arrives at 



W{Ai,A2) = exp 



A. G' 



Y,Pr{A2,A) exp[A2 G\zi) 



(44) 



where /3i's are sums of polynomials in A2 times functions of A. One can write 
the above equation in a form more manifestly-symmetric with respect to Ai and 
A2: 



W{Ai,A2) =exp 











exp 




+ exp 





Y,p,{A2,A) exp[A2G'(z,)], 



(45) 



where 



A -exp [AG' (-i)] , 



(46) 



and use has been made of the fact that G' is an odd function. 

For YM2 , one can perform a more explicit calculation. Here it is known that 
for A> Ac^ the density is equal to one for — & < z < 6, where h is some positive 
number less than a. Also, it is seen that a and 6 both tend to ^ as A tends to 
00. So there are only two z^'s and one has 



H{z) 



lni±| 



a{A) 



1 



1 



(z+i)2 



+ 0(«'), 



(47) 



where 



1 /I 



2 V2 



dy y 



p{v)- 



(48) 



(One can see that H is an odd function and its residue at 2 = should be equal 
to one, so there are no first-order poles at z = ±5-) Putting this in (|13|l . one 
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arrives at 



W = - 



dz z 



"2,111 z + \ 



1 — a 



(^_l)2 (,+ 1)2 



= expl-- 



1 + a ( -1 - ^2 + ^ 



3 

2 



A, 



+ a exp ( ^ ) + 0(a2). (49) 



It is known that for A2 = A, Wg should become unity. This shows that up 
to the leading order, 

A' 



a = exp ( - — 



(50) 



So, 

W = exp 



A2 
' 2 



1 + exp 



-1-^2 + ^ 



+ exp 



' 2 



+ 



This is the same relation obtained in [24], using a different method. 



(51) 
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